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A SENSITIVITY ANALYSIS OF VOLCANIC AEROSOL DISPERSION 


IN THE STRATOSPHERE 
By 

Carolyn F. Butler 
ABSTRACT 

A computer sensitivity analysis has been performed to determine the 
uncertainties involved in the calculation of Volcanic aerosol dispersion in 
the stratosphere using a 2-dimensional model. The Fuego volcanic event of 
1974 was chosen for this study. Aerosol dispersion processes that were in- 
cluded are: transport, sedimentation, gas phase sulfur chemistry, and 

aerosol growth. Calculated uncertainties are established from variations in 
the stratospheric aerosol layer decay times at 37N latitude for each dis- 
persion process. Model profiles are also compared with lidar measurements 
obtained at the NASA Langley Research Center in Hampton, Virginia (37N). 

Results of the computer study are quite sensitive (factor of 2) to the 
assumed volcanic aerosol source function and the large variations in the 
parameterized transport between 15 and 20 km at subtropical latitudes. 
Sedimentation effects are uncertain by up to a factor of 1.5 because of 
the lack of aerosol size distribution data. Finally the aerosol chemistry 
and growth, assuming that the stated mechanisms are correct, are essentially 
complete in several months after the eruption and cannot explain the differ- 
ences between measured and modeled results. 
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1., INTRODUCTION 



A characteristic of the earth’s stratosphere Is a dust layer attributed 

to volcanic eruptions. This dust layer is made up of liquid or solid aerosols 

and normally has its peak particle mixing ratio between 18 and 21 kilometers 

(11) . Stratospheric aerosols are predominately sulfate particles"— possibly 

sulfur dioxide converted through a series of reactions to sulfate and then 

hydrolyzed to sulfuric acid (5). While the major portion of all stratospheric 

aerosols is due to volcanic eruptions, other natural and anthropogenic aerosols 

may contribute to the layer via troposphere-stratosphere exchange processes 

(24). During periods of low volcanic activity, background concentrations of 

3 

aerosols are on the order of 0.5 particles per cm or less with sizes ranging 
from about 0.1 to 1.0 micrometer in radius (12). Particles larger than 1.0 
micrometer fall out rapidly by sedimentation and particles smaller than 0. 1 
micrometer grow by condensation and coagulation processes (30). 

While pollutants in the troposphere are quickly dispersed by the winds or 
removed by rainout (or washout) , stratospheric transport processes are weak 
and stratospheric gases and particulate matter have much longer residence 
times. Furthermore, this region is close to radiative energy balance and 
introduction of foreign gases or particles could disrupt the radiative energy 
budget resulting in climate modification at the earth's surface. The aerosol 
layer is of particular importance because it resides in the stratosphere. 
Aerosols affect radiation by both their absorption and scattering properties 
which are determined, in turn, by various aerosol characteristics such as 
size, shape and composition. The properties of the stratospheric aerosols 
and their effect on the radiation balance have been reviewed by Cadle and 
Grams (1) . 



In order to adequately assess possible climatic Impacts , more must be 
known about the dispersion and residence times of aerosols following volcanic 
eruptions. This work considers the dispersion of aerosols from one particular 
event — the eruption of Mfc. Fuego In Guatemala (15N) on 14 and 17 October 1974. 
Atmospheric transport, particle sedimentation and aerosol physical chemical 
processes are evaluated for that event by comparing results with observations. 

Remote sensing techniques have been used to monitor stratospheric aerosols 
from ground stations and from aircraft. One such technique, the laser radar 
(Ildar), has been successfully utilized since 1963 to define vertical profiles 
of aerosol layers. Briefly, the lidar technique consists of a laser transmitter 
which emits a pulse of light vertically into the atmosphere r where the incident 
photons are absorbed and scattered by molecules and aerosols. The 180° backscat- 
tered light from both the molecular atmosphere and the aerosols is collected by 
a telescopic receiver located co-linearly with the transmitter, The principle 
of the lidar and the lidar calibration is described in more detail by Northam, 
et al. (21) and Remsberg, et al. (33). 

A lidar measure of the aerosol mixing ratio is the scattering ratio R 

b 

R » 1 + f /f (1) 

s a m 

where f and f are the aerosol and molecular backseat tering functions, respec- 
a m 

tively (21). The f-values are products of the species cross-section and number 
density. The portion of any R -value greater than one represents backscattering 
from aerosols. If the aerosol cross-section is constant with height, the scat- 
tering ratio profile is a direct measure of aerosol number density. For an aerosol 
size distribution which remains constant with time, the stratospheric aerosol 



scattering ratio can be used to access the relative change in aerosol number 
density. 

Figure 1 shows two R profiles obtained by lidar at Hampton , Virginia. 

S 

The plot for January 2, 1975, represents an enhanced aerosol layer due to 
volcanic activity in Guatemala in October 1974 (25). The February 19, 1976, 
profile resembles a near- background aerosol level and shows the depletion of 
the January 2 layer over 13 months. As a result of the substantial quantity 
of lida.. 1 data available from Hampton, Virginia (31) and from facilities at 
other latitudes, an analysis of theoretical models of the latitudinal and 
vertical dispersion of the stratospheric aerosol layer can be conducted. 

The present study involves a sensitivity analysis of the aerosol layer to 
the various assumptions which are, made for the aerosol model in the disperison 
calculations. Section 2 describes the circulation model and the computational 
procedure. The aerosol layer processes and the assumptions concerning the 
aerosol model are explained and justified in Section 3. Gas phase chemistry 
and aerosol growth effects are also discussed there. Section 4 then describes 
results of the aerosol sensitivity studies and the validity of those results. 
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Figure 1. Aerosol scattering ratio profiles from lidar at Hampton 
Virginia, following the eruption of Volcan de Fuego in 
October 1974. 


2. MODEL DESCRIPTION 


The aerosol distribution is predicted in the form 

N(t) - N(t 0 ) + (dN/dt) At ^ 

where N(t Q ) represents the Initial aerosol distribution (mass mixing ratio) 
and At is the time step. The rate of change of aerosol mixing ratio dN/dt 
is defined by 



where (SN/3 1) tr Is the transport terra, t> Bed Is the sedimentation term, 

(8N/9t) Is the aerosol growth rate and the gas phase chemistry term. 

8 r 

For simulations presented in this paper, the transport is specif if led by 
monthly mean winds and eddy diffusion parameters derived from the seasonal 
circulation in Louis' Model II (18, 2). The model extends from 0 to 50 km 
in altitude with a vertlcle grid spacing of ,1 1cm and from 90N to 90S with a 
horizontal grid spacing of 5 degrees. The continuity equation for the aerosol 
mass mixing ratio is integrated at specified time steps using a semi- imp licit, 
centered-dlf ference scheme. Louis' model has approximated reasonably the dis- 
tributions of trace gases and radioactive debris in the stratosphere. In 
particular, the analysis of a volcanic event represents dispersion from a point 
source, similar to that for radioactive bomb debris. 

The latitudinal boundary conditions are imposed such that there is no flux 
at the poles. At 50 1cm the boundary condition allows mass to be advected out of 
the model, but not diffused out, however, this boundary is well above the aerosol 
layer and should have little effect on the aerosol dispersion. Ac the lower 
boundary, one has the choice of either specifying a constant mass mixing ratio at 
the boundary or a constant flux through the bouvidary. 
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Thn sedimentation term Is simulated by applying fall speeds for various 

-3 

aerosol sizes and a density of 1.5 g cm . These aerosol fall speeds have 
been tabulated by Kasten (14) for several particle sizes and at various 
altitudes. The assumptions and techniques. used to predict the effects of 
sedimentation will be discussed In detail in section 3. 

The gas phase chemistry Is considered In terms of its Importance to 
the aerosol mass loading, and subsequent aerosol growth has been incorporated 
Into the model as >.11. The effects of the chemical and growth term will be 
discussed in sections 3 and 4. 

To initialize and run the model, one must specify an initial gas and 

aerosol mass distribution, a background gas and aerosol mass distribution,. 

the time step and the lower boundary condition. Initial conditions of aerosol 

mass distribution N(t ) for the model are estimated from several data sources. 

o 

HiJijvi resolution infrared and visible satellite photographs from SMS (released 
by the National Oceanic and Atmospheric Administration in Rockville, Maryland) 
and from the DMSP (U. S. Air Force satellite) are employed to estimate the 
eruptions of October 14 and 17, 1974 (32). Dust clouds were observed to move 
to the northeast and to the west. An analysis of rawinsonde profiles obtained 
at Swan Island (off the east coast of Guatemala) shows winds from the southwest 
up to about 15 km with a sharp reversal to winds from the east above 16.5 km. 
Trajectory analysis of dust layers in these two flow regimes explain some of the 
subsequent observations at higher latitudes (25). 

An analysis of the photographs and the local wind profiles indicate that 
the October 17 event was responsible for the bulk of the 20 km layer; a grid 
based on that event is used in the model. The model is started on October 19 

b 



at 1200Z to allow for some spread of the cloud. Amounts of Injected ash and 
sulfuric acid gas are estimated as a factor of 20 less than the Initial distri- 
bution for Agung as estimated by Cadle et al. (2). (Cadle, et al. determined 
the Fuego e\*ynt to be only a factor of five less than Agung). Data taken by 
iidar at Hawaii (4) on October 29 are applied to verify the Initial vertical 
profile of the dust layer. The vertical width of the layer at half maximum as 
computed by the model for October 29 at 20N is compared with the observations 
(half-width of 0.8 kra) at Mauna Loa Observatory in Hawaii for that dat^. The 
shape of the model profiles at all altitudes where aerosol had been transported 
in that. 10-day period is adjusted to agree with a zonal average of the Hawaii 
observations. These adjusted profiles then represent the initial conditions on 
October 29 for the aerosol source and are shown in figure 2 . The sensitivity 
of the calculations to the specification of the source function is discussed in 
a subsequent section. 

The choice of the initial aerosol size distribution is important for the 

aerosol sedimentation and growth processes and that distribution is shown in 

Figure 3. The rationale for this distribution is given in Appendix I. The 

-10 

background aerosol mass mixing ratio is estimated to be 2 x 10 and the 
background sulfuric acid gas concentration is estimated from Figure II-4 in 
Appendix II. 

A time series of dustsonde profiles after the eruption of Fuego shows 
considerable variation in the adrosox mass mixing ratio in the upper tropos- 
phere over Laramie, Wyoming (8). This indicates that mixing of aerosol from 
the stratosphere to lower altitudes can occur in a non-uniform manner. As a 
result, the entire tropospheric and stratospheric circulation of Louis' Model 
II is used and a constant mass mixing ratio of 2 X 10 ^ is assumed for the 


7 





PERCE 



0.1 0.3 0.5 0.7 0.9 

RADIUS (micrometers) 


"Figure 3. Initial size distribution estimate of 
Fuego event used in this model. 
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model at a lower boundary of 2 kilometers. 


This time step is chosen such that it is at least an order of magnitude 
less than the transport relaxation time. This insures stability in the mathe- 
matical computations of the continuity equation (eq» 3) (18). The transport 
relaxation times for the mean winds and eddy diffusion parameters are defined 
as follows 


t « H /v; t « H /w; t v * H /K . 
v y w z' K y yy; 


» H /K ; -HH K 

K. z zz K z y ; 

zz yz 


yz 


where v and w are the latitudinal and veritcal wind components, respectively; 
y and z refer to the latitudinal and vertical directions, respectively; the 
K-values are the eddy diffusion coefficients, and; H is the aerosol scale ( 
distance over which the mixing ratio changes by l/e. Initially, the aerosol 
scale height is small due to the steep graidients of the source, so one hour 
time steps are used. As the dust disperses, the time step is gradually in- 
creased to 12 hours. 


* 
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H. AEROSOL LAYER PROCESSES 


This section deals with each of the processes in equation 3. Even though 
the gas phase chemistry Is not included in the model calculations, assumptions 
concerning this process are discussed and their effect on the aerosol dispersion 
is estimated. 

A. Transport 

The basic equation for time dependent tracer transport is 

B. m _ v • N + ~ V • <pK • VN) (4) 

at p 

where p is density of air, V is velocity, V is a gradient operator, K is a second 

order diffusion coefficient tensor, and N is concentration of the tracer species. 

♦ 

The first term on the right of (4) represents advection by the mean winds while 
the second term represents gradient diffusion. Source and sink terms are expressed 
separately in equation (3) . 

The transport term in this modeling attempt is represented by monthly mean 
winds and eddy diffusion parameters derived from Louis' circulation Model II 
(18, 2), which adopt the mean winds derived from observations by Newell, et al. 

(20) below 15 kilometers. Above 15 km the mean winds are computed using the 
thermodynamic and continuity equations. The eddy diffusion parameters are re- 
presented by a tensor coefficient due to physical mixing in the atmosphere and 
also to the mathematical terms when zonal averaging is performed. The observed 
distribution of ozone and its observed flux at the ground were used by Louis to 
adjust the diffusion coefficients with the mean circulation. 

The initial assumption in this report is that Louis' Model II has reasonably 
predicted the dispersion of radioactive bomb debris and therefore should also ap- 
proximate the dispersion of the volcanic dust layer. Even though this particular 


11 



model possesses some deficiencies (largely due to a lack of data sets for para- 
meterizing the circulation), It should be adequate for sensitivity studies of 
the source function, and of the sedimentation, chemistry and growth terms of 
equation (3). The sensitivity of the results to the uncertainty in the trans- 
port term Itself is not addressed explicitly in this study. 

It should be noted that the choice of the boundary condition can change 
the net effect of the mean winds and large-scale eddies. The sensitivity of 
the model to the boundary condition is tested by varying the boundary from 2 
kilometers to 10 kilometers. These results will be discussed in Section 4. 

B. Sedimentation 

Hunten (9) has discussed the importance of aerosol sedimentation rates for 

determining the residence*' times of volcanic aerosol layers. Aerosol fall speeds 

for Various aerosol sizes tabulated by Kasten (14) for particle densities of 1.5 
-3 

g cm are applied to the initial aerosol size distribution derived in Appendix 
I and plotted in Figure 3. This is a best estimate of the size distribution. 

An upper limit on the sedimentation rate is obtained with the April 11, 1963 size 
distribution measurements by Mossop (19) after the eruption of Mt. Agung. Mossop's 
distributions overestimate the number of larger sized particles from the Fuego 
eruption since the Fuego event was not as violent as the Mt. Agung eruption. 

Mossop* s impactor was also biased against the smaller particles. Nevertheless, 
this data set should at least give an upper limit to the sedimentation rate. A 
lower limit sensitivity to the sedimentation rate is also tested by applying ap- 
propriate fall speeds to a log-normal size distribution representative of a back- 
ground size distribution of aerosols (23). 
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The total initial size distribution is divided into five size ranges with 
mean radii of 0.10, 0.28, 0.50, 0.69 and 0.90 micrometers. Distribution of mass 
mixing ratios of five different sizes are then calculated as functions of alti- 
tude and time. 

C. Gas Phase Chemistry 

The gas phase chemistry of the SO^ to aerosol conversion is initially con- 
sidered with a one-dimensional model to determine its importance in the long-term 
aerosol dispersion. A simple S0 2 to H 2 S0^ gas phase chemistry is assumed to follow 
the route 

so 2 -*■ hso 3 + so 3 •* h 2 so 4 . 

The specific reactions considered are listed in Table II-I of Appendix II. Junge 
(13) lists several other reactions that are likely to occur but since the rates of 
these reactions are unknown * 3 nd the reaction schemes have not been confirmed, they 
have not been included. A steady-state model of aerosol gas phase chemistry via 
(5) is described in Appendix II. The S0 2 -> conversion time can be estimated 

by examining the photochemical relaxation times of reactions 1 through 4 of Table 
II-I. The photochemical relaxation time is defined as follows 

*,h - 1/K * <6) 

where K is the effective reaction coefficient. The effective rate is the product 
of the number density of the reactant species and the reaction rate. Several as- 
sumptions must be made to calculate the effective rates. The reaction rates are 
not always well known and can vary by orders of magnitude from one reference to 
another. The concentrations of the reacting gases are sometimes uncertain as well, 
so one must assume a particular concentration profile. The calculations of the 
photochemical relaxation times are further complicated by variations of the reaction 
rates with altitude and time of day. Reaction rates will vary with changing 
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temperatures and/or the presence of sunlight and certain species concentrations 
can change by orders of magnitude as a function of altitude. The rate coeffi- 
cients assumed in this work are listed in Table II-I with references. Park and 
London (22) have developed a photochemical model including profiles of the species 
OH, M, 0, and H0 2> Their species number densities at 20 km are used to estimate 
the photochemical relaxation times. 

Each phase along the S0 2 ->■ HSO^ ■> S0 3 •* H 2 S0^ route appears to be quite rapid. 
The rate determining reactions in the S0 2 phase (see Appendix II, Figure II-I) in- 
dicate that photochemical relaxation times are on the order of 1 to 10 days. If 
the uncertainties in the reaction rates and species concentrations are taken into 
consideration, the relaxation time could be as much as 50 days. 

If the amount of S0 2 injected by the Fuego volcano into the stratosphere far 
exceeded the amount of aerosol injected, then the chemistry term would contribute 
a major bulk of the aerosol mass to the ambient background layer. However, judg- 
ing from estimates of Gadle, et al. (2),. the eruption cloud contains about equal 
amounts of SC> 2 (or gaseous sulfuric acid) and particulate matter. Since the gas 
phase chemistry appears to be very rapid, its effects should be small and the 
chemistry has not been included in the long term dispersion model for calculations 
of aerosol decay over the order of a year. 

D. Gas to Aerosol Conversion and Growth Effects 

Aerosol growth by both coagulation and net condensation mechanisms has been 
evaluated by Turco, et al. (Id) and Hamill et al. (30). Coagulation processes 
cause a change in the size distribution with time, where larger particles grow 
at the expense of smaller ones. Hamill, et al. (30) show that coagulation is 
not important for particles with radii greater than 0.2 micrometer. Therefore, 
coagulation, a process which conserves aerosol mass, is not included in these 
model calculations. 
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Condensation, or sulfuric acid gas to aerosol conversion would actually 
add to the total aerosol mass. The condensation process varies directly with 
sulfuric acid gas number density and, for a volcanic event, would be a function 
of altitude. For a layer of aerosol and gas injected at 20 km, the net effect 
of growth by condensation on the aerosol profile would be to prolong the exist- 
ence of the aerosol layer. 

To simulate the condensation process the growth rates calculated by Hamill 
(30) are used. These rates are given as a function of altitude and particle 
radius and they are normalized by the number of molecules of gaseous sulfuric 
acid present. The parameterization of the growth rates is simplified by choos- 
ing those values at 20 km only. Two different sulfuric acid gas profiles are 
tested with the model— gas mixing ratio values at the source equal to and twice 
the magnitude of the aerosol source profile (based on estimates by Cadle, et al. 
(2) of about equal amounts of ash and sulfuric acid injected by Agung) . At grid 

points other than where the volcanic source occurred, a background mass mixing 
-15 

ratio of 6 x 10 is used (see Appendix II, Figure II-4). At each time step and 
grid point and for each of the Initial five mean particle radii, particle growth 
rates are calculated. The change in the mass due to growth is added to the aero- 
sol mass and simultaneously subtracted from the sulfuric acid vapor mass. The 

sulfuric acid vapor profile is not allowed to fall below the equilibrium value of 
-15 

6 x 10 since it is assumed that the chemical processes (not included in the 2- 
dimensional model) maintain at least some background amount. The particle size 
is then updated and a corresponding sedimentation rate is calculated. Results 
of the growth sensitivity studies are presented in Section 4. 

E. Relationship Between Calculated and Observed Mass Densities 

The model aerosol distribution is calculated in terms of mass mixing ratio 
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and this quantity is compared later with lidar backscatter measurements froti the 
aerosol layers. The lidar backscatter results from the product of aerosol number 
density (which is proportional to aerosol mass) and aerosol backscatter cross 
section. If the backscatter cross section is constant with time, the trends in 
the lidar backscatter can be directly compared with the model results. 

The aerosol cross section depends on the refractive index, particle shape 
and size distribution. Measured aerosol samples are predominantly pure liquid 
spheres (probably sulfuric acid) or ash nuclei coated with liquid, and since the 
refractive index is determined principally by the outer portion of the aerosol 
particle, no change in refractive index with time is assumed. Measurements by 
Hofmann and Rosen (8), however, indicate a changing aerosol size distribution, 
at least in the early stages of the volcanic dispersion cloud. Remsberg, et al. 
(33) have analyzed the sensitivity of ruby lidar backscatter returns to various 
size distribution models and they have determined that this effect is only im- 
portant for background aerosol size distributions, not volcanic distributions. 
Therefore, direct comparisons between aerosol mass mixing ratios and lidar back- 
scatter should be valid at least for time periods up to one year after the Fuego 
eruption. 
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4. RESULTS OF SENSITIVITY ANALYSES 


A. Transport 

The sensitivity of 9N/3t to various transport models has not been tested 
here explicitly even though transport is a dominant term in equation (3) . This 
analysis might be achieved by inputting the Fuego dust source into Louis' other 
model, Model I (18), in which the mean winds are greater by a factor of 2, or 
into other available circulation models. Glatt and Wldhopf (40), in particular, 
have applied different radiation and turbulent diffusion models to the calculation 
of the seasonal meridional circulation patterns according to the formulation of 
Louis and have tnen compared their results with those of Louis. That comparison 
yielded large differences in the tropical and subtropical lower stratosphere, 
just where the Fuego dust was Inserted. In general, their new circulation fields 
indicated reduced vertical transport at 20 km and a delay in the meridional trans- 
port to southern latitudes. It is also possible that the 1974-75 seasonal circu- 
lation patterns exhibit differences from those developed from earlier atmospheric 
data sets by either Louis or Glatt and Widhopf. Due to these differences, circu- 
lation patterns and observations of past volcanic eruption clouds are briefly 
examined instead to determine whether Che transport model compares with the general 
predictions of the Fuego dust route. 

According to Lamb (15), the prevailing zonal winds can quickly carry the dust 
layer around the globe. A typical circuit would take from two to six weeks depend- 
ing on the latitude of the source. The zonal components of these winds axe several K 

orders of magnitude greater than the meridional and vertical components. Therefore, 
spread of the dust to other latitudes and altitudes should be much lower. The non- 

l 

symmetric wave structures of these zonal wltids, on the other hand, transport the 
dust particles in the north-south direction. These processes are represented by 
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the mean meridional circulation and the large scale eddy diffusion as discussed 
earlier. The three cell meridional circulation Is clearly visible In Louis' 
model of the winter circulation (Figure 4) . Eddy transfers which operate In 
both directions, appear to dominate the mean circulation In raid-latitudinal 
spread In winter (see Park and London (22)). 

Observations Indicate that in late autumra there Is a sudden spread of dust 
Into higher latitudinal belts. Dust from the Krakatoa (Indonesia) (May and 
August, 1883) and Agung (Ball) (February and March 1963) eruptions, both In 
equatorial latitudes, spreads quickly to about 35N and 35S but was not observed 
at higher latitudes until late in the following autumn of each hemisphere. In 
fact, there Is evidence that tnis spread of volcanic dust Into other latitude 
zones is made during the great seasonal circulation changes. 

Since the Fuego event occurred In the Northern Hemisphere during the autumn 
season, one would expect the dust to arrive at the latitude of Hampton, Vir- 
ginia (37N) in just a few weeks. Observations of past volcanic events shows 
this to be the case. The eruption of Mt. Agung (8S) occurred during the autumn 
season of the Southern Hemisphere. The Agung dust was first detected over Mel- 
bourne (38S) some 30° to the south In five to seven weeks. Maximum concentrations 
were observed after about 4 to 6 months at Melbourne but were not observed until 
one year later after the- eruption at latitudes 40-45S where reverse circulation 
develops. The Fuego dust had only to travel 22 degrees to reach Hampton and, In 
fact, initial sightings were made in Hampton in 4-6 weeks and maximum concentrations 
were observed some three months after the eruptions. If we allow for the greater 
transport distance from Agung to Melbourne compared to Fuego to Hampton and for 
the lower latitude of the source (8S vs 15N) , then the initial appearance of the 
dust and the time of maximum concentrations compare very well. Remsberg and 


18 



ALTITUDE, KILOMETERS 



90 60 30 0 30 60 90 

NORTH SOUTH 

DECEMBER - FEBRUARY 


12 -1 

Figure 4. Integrated mass flux (10 g sec ) from 
the winter mean circulation of Louis ' 

Model I. The dot represents the altitude 
and latitude of the initial Fuego volcanic 
cloud. 
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Northern (25) have used lower stratosphere circulation maps for October 1974 to 
explain In more detail the latitudinal spread of the Fuego dust layer. 

Figure 5 displays the Integrated aerosol mass density from Ildar observations 
at Hampton, Virginia (37N) between 16 rnd 21 km as a function of time after the 
eruption. Considerable variability Is present In the early returns. The vari- 
ability in the Ildar returns through mid-December represents longitudinal inhomo- 
geneltles of the volcanic dust, a feature that cannot be simulated by a zonally 
averaged model* The Integrated Ildar data are obtained by summations over 1 km 
altitude (z) Increments of 


(R s (z) - 1) N m (z) (7) 

where R (z) Is the average scattering ratio and N (z) Is the molecular number 
s m 

density (see eq. (1)), * Thus, the relative aerosol column density applies to a 
5-km column of 1-cm 2 cross section. The quantity on the ordinate is then equi- 
valent to (a h ) N (z) where a and a are aerosol and molecular cross sections, 
respectively, and (z) is the aerosol number density which is directly propor- 
tional to aerosol mass density, M (z) . Therefore, the dimensions of the ordinate 
in Figure 5 are ML'* 2 . No attempt has bean made to actually compute absolute mass 
densities from the Ildar data. That is, no adjustment has been made for possible 

variations of a with time or altitude, but theoretical studies by Remsbcrg, et al. 
a 

(33) indicate no appreciable changes in c & for the ruby lidar data. 

Figure 6 displays the model integrated aerosol mass density between 16 and 
21 km for an average of 35N and 40N latitude. Only the effects of the transport 
term are considered here. Two curves are generated to determine the sensitivity 
of the model calculations to the lower boundary condition. A constant aerosol 
mass mixing ratio of 2 x 10 _1 ° was imposed at 10 km and then on a second run at 
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10 km. J3y setting the boundary by the lower altitude, a decrease occurs in the 
vertical concentration gradient and the flux across the boundary and into the 
troposphere, thereby maintaining more mass in the 16 to 21 Urn column. A 23 
percent difference is noted between the two runs one year after the Fuego erup- 
tion. The preferred curve is the one for the 2 km boundary. The justification 
for this choice is based on dustsonda measurements by Hoffman and Rosen (8) 
indicating a variable rather than a fixed mixing ratio at 10 km altitude from 
three months after the eruption. This boundary condition also requires a 
coupled stratosphere- troposphere model such as that which has been applied in 
this work. 

The time of the occurrence of the measured maximum aerosol load at 35N is 
simulated very well by the model. If the magnitude of the lidar data in Figure 
5 and the upper model curve in Figure 6 are normalised to their peak values, the 
respective time rates of change of the aerosol column are 8 months and 4.5 months. 

Figure 7a displays aerosol profiles for 35 to 40N for February and May 1975 ; 
the Ildar data are plotted in terms of aerosol scattering ratios (a pseudo-aerosol 
mixing ratio) while the dustonde profiles from the University of Wyoming (8) are 
in terms of aerosol number density mixing ratios for particles greater than 0.15 
micrometer in radius. Although there are some amplitude variations between the 
Ildar and dustsonde data, the mean altitudes of the layer peaks and the widths at 
half-maximum are comparable. The corresponding model profiles due to only the 
transport term ore presented in Figure 7b and the model profile wi' ; th at half- 
maximum is overestimated, which means that either the transport is too rapid or 
that the source function has been specified Improperly. This overestimation is 
also evident in comparison between model profiles and lidar data at other latitudes. 
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Figure 7a and 7b. Vertical profiles of aerosols from observations 
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If the vertical spread of the model aerosol layer is too rapid then aerosol 
material will be transported both above and below the 16 to 21 km column. The 
Ildar data do not show significant aerosol above 21 km. The effect of enhanced 
vertical transport of aerosols by the model Is presented In Figure 8 by plotting 
curves for 16 to 21 km and for 16 to 25 km. One year after Fuego, about 25 per- 
cent of the aerosol between 16 and 25 km is above 21 km and the 1/e decay rate 
is slightly different for the two curves. However, inclusion of the sedimentation 
term effectively eliminates the spread of material above 21 km and therefore the 
16 to 21 km column density curve is selected for all further sensitivity studies 
in this work. 

In conclusion, the transport term appears to represent the meridional trans- 
port between 15N and 37N quite well, while overestimating the vertical transport. 
The uncertainty in the vertical transport which may be due partly to the speci- 
fication of the source function in the model, has a significant effect on the 
aerosol dispersion calculation. By inserting the same aerosol source into another 
2-dimensional model with the same grid and time step, one could conduct a further 
check on the uncertainty due to the parameterized circulation as well. That effort 
is outside the scope of this study. 

B. Source Function 

The proper specification of a point source event in a zonally averaged model 
that incorporates gradient diffusion is extremely important but difficult. Louis 
(18) alluded to this problem when he was formulating his 2-dimensional model and 
he was forced to consider observations of ozone in order to improve the parameteri- 
zation of his model. Cadle et al. (2) then applied the model to another point 
event in the lower stratosphere, the Agung eruption. They adjusted their source 
function so that the model results agreed with subsequent observations of the 
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dust layer at 30N and 40N, latitudes which are far to the north of the source 
(8S). Their approach required that they assume the circulation to be correct, 
because the circulation and source function affects could not be decoupled. 

Thus it is no surprise that their modeled aerosol decay curves agreed with the 
mid- latitude lidar measurements. Then they further assumed for the Fuego event 
that the source function was the same as that for Agung, except that the source 
latitude was moved from 8S to 15N. 

To test the sensitivity of the model results to the source function in the 
present study, the circulation has been assumed correct and several source func- 
tions have been tried. The source grid in Figure 2 is thought to have the correct 
shape but the observations against which it was checked are taken from only one 
location, Hawaii. Therefore, the zonal averaging process can impart errors to 
the profile shape, and more Important, to the magnitude of the peak mixing ratio, 

The shape of the initial vertical profile has been varied, while retaining 
the same peak mixing ratio, but there are no appreciable differences in the result- 
ing decay curves. No runs were made to test the effects of significantly different 
source function magnitudes. Although the vertical profile is not resolved very 
well by the 1 km model grid (39), the dispersion appears to be numerically stable 
because of the inclusion of the diffusive transport term in the model. 

The sensitivity of the advective and diffusive dispersion rate to the assumed 
initial latitudinal source gradient (see equation 4) has not been tested here, 
but from Louis’ (18) studies it is one of the important uncertainties t The lati- 
tudinal gradient from the model run for 60 days after the start date is compared 
with that published by Cadle et al. (2) and their concentration gradient is 
weaker. This fact may partially explain why their model decay times are longer 
than those in the present study. Because the details of their source are not 
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known, further comparisons with their results are not warranted. In conclusion, 
the uncertainty in the source function for a 2-dimensional model can only be 
minimized by waiting about three months after the eruption before initiating 
the model at which time the longitudinal Inhomogeneities are smoothed out and 
measurements of the latitudinal and vertical concentration gradients are available. 

C. Sedimentation 

Figure 9 represents the integrated aerosol density between 16 and 21 km as 
a function of time after the eruption for the case where the sedimentation term 
is included. The solid curve represents the case where fall speeds were applied 
to a log-normal background type aerosol distribution (23) of mean radius 0.0726 
micrometers. The implication is that the sedimentation term has no effect on a 
background- type size distribution. The triangles are results for the April 11, 

1963 volcanic aerosol size distribution taken from Mossop (19). This distribution 
contains many large particles which settle out quickly and as such represents an 
upper limit on the sedimentation effect. The circles represent results for the 
"best estimate" Fuego initial aerosol distribution shown previously in Figure 3. 
This distribution is used in the aerosol growth sensitivity studies in the rext 
subsection* 

The. principal effect of the sedimentation term is to decrease the amplitude 
of the total curve to shorten the decay time to 2.5 months, and to shift the peak 
value to the left by about ten days. By not including this term (compare solid 
curve with circles), one would overestimate the aerosol mass and calculate incor- 
rect decay times. The triangles approach a lower limit in Figure 9, because of 
the 2 s 10 ^ lower limit in mass mixing ratio. 

D. Aerosol Growth 

Figure 10 presents relative aerosol mass density curves for two conditions 
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1) no gas to aerosol conversion and 2) condensational growth assuming an 
initial sulfuric acid vapor profile that is equal in magnitude to the sulfuric 
acid in aerosol phase. Upper limit estimates of the ratio of vapor to aerosol 
phase sulfuric acid were obtained from Cadle et al. (2) and condensational growth 
rates were taken from llamill, et al. (30), Although a growth rate at 20 km was 
assumed, the rates at the other altitudes between 15 to 25 km vary from it by less 
than 15 percent for particles with radii between 0.1 and 1.0 micrometers. For a 
fixed altitude the additional uncertainties in the growth rate equation derived 
from Hamlll et al. (30) have not been modeled, but they are believed to be small 
because Hamill et al. were able to characterize the background stratospheric 
aerosol layer quite well using the same expressions. The dashed curve in Figure 
10 is a "best estimate" of the effects of growth on the volcanic decay curve and 
is based on initial sulfuric acid gas to aerosol ratio of one- third. The basis of 
that estimate is discussed later in this section. 

Figure 11 presents the changes in the initial aerosol size distribution curve 
(Figure 3) at 20 km for three latitudes — 35N, 15N, and 25S. Transport, sedimenta- 
tion, and growth are combined to give these curves. At 15 N one finds that the 
effects of aerosol growth are essentially completed after one month. At 35N the 
effects of sedimentation compensate those from aerosol growth during the early 
months and a loss of larger particles by sedimentation is apparent one year later. 
At 25S only the sedimentation process effects the size distribution, because the 
volcanic gaseous sulfuric acid is essentially depleted before reaching this 
latitude. The log-normal background distribution is plotted on that same figure 
for comparison. 

Table 1 lists the number and mass percentages for each size range for both 
the log-normal and the initial volcanic distribution. The bulk of the aerosol 
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Figure 11. Changes in the initial size distribution curve 
at 20 km for three model latitudes (35N, 15N, 
and 25S) . 
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mass occurs at the larger sizes for a volcanic situation. Therefore, any error 
due to the collection efficiency of the aerosol Instruments at 0.1 micrometer Is 
minimized when calculating total aerosol mass density. As a check to the computed 
effects of condexisational growth on the aerosol size distribution curve, Hamlll 
et al. (30) estimate a shift of the curve to larger particles by 0.1 micrometer 
In about 6 months. Figure 11 (where both condensation and sedimentation are act- 
ing) shows that at the latitude of the source, a shift to the right of 0.14 
micrometer occurs In one month with little change thereafter. At 35N a shift of 
0.07 micrometer Is apparent after one year. 


Table 1 

Stratospheric Aerosol Size and Mass Distribution 


Initial Volcanic Distribution 

Log-Normal Background Distribution 

Radius 

Percent by 
Number 

Percent by 
Mass 

Radius 

Percent by 
Number 

Percent by 
Mass 

0.10 

24 

0.4 

0.10 

94 

30.2 

0.28 

52 

19.7 

0.28 

5.4 

38.1 

0.50 

18 

38.9 

0.50 

0.5 

20.1 

0.69 

5 

28.4 

0.69 

0.09 

9.6 

0.90 

1 

12.6 

0.90 

0.009 

2.1 


A comparison has been made in Figure 11c between the model results for Fuago 
and the observations by Mossop for Agung, both curves for one year after the res- 
pective eruptions. Note that the size distribution slope between 0.1 and 0.28 
micrometers is greater for the Agung data. This difference exists even though one 
might expect the opposite because it is probable that there were proportionately 
more large particles associated with the more intense Agung event. Also by con- 
sidering the reduced collection efficiency at 0.1 micrometer for Mossop' s impactor, 
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the slope of the observed curve would be even steeper than shown in Figure 11c. 

To explain this discrepancy, more 0.1 micrometer radius particles must be gene- 
rated in the model. The calculations presented here do not allow formation of 
new particles by nucleatlon and their subsequent growth by coagulation and 
condensation to the 0.1 micrometer size range. Nucleatlon, then, represents a 
possible important process that could be Included in the present model if the 
rate of replenishment of 0.1 micrometer particles could be estimated, Inclusion 
of this process, however, should have little effect on the resultant aerosol 
decay curve unless the formation of sulfuric acid vapor is much slower than es- 
timated in Appendix II. A test has been performed where new particles (0.1 
micrometers in radius) are allowed to form at a rate equal to the loss due to 
growth of the particles in the 0.1 micrometer range on the assumption that these 
smaller particles have a self-preserving property. These "nucleated'’ particles 
are then dispersed by the model transport processes but they are not allowed to 
grow by condensation. As expected, there is no effect on the aerosol decay curve. 
The size distributions in Figure 11 change considerably, however, with the largest 
number of particles now in the 0.1 micrometer range at all latitudes. 

Model profiles for February and May 1975 at 37N are presented in Figure 12 
for the sedimentation and for the sedimentation plus growth cases. These model 
profiles should be compared with the measured and modeled profiles in Figure 7. 
Note that the altitude of the maximum aerosol mass mixing ratio is below the 16-21 
km range for both cases. A large percentage of the total stratospheric aerosol 
mass has fallen below the 16 to level and, as a result, the model layer decay 
time is about 2.5 months. The measured profiles in Figure 7a do not exhibit this 
rapid depletion of aerosol. The principal reason for this discrepancy is beleived 
to be the transport in the lower stratosphere in the model. 
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The rapid growth process , which wan complete after several months, agrees 
with the, qualitative, conclusions from s tud 1. •, by Hirono et al. (41) and Hoffman 
and Rosen (8), and from the sulfur isotope, data of Gastlemau (42) for another 
eruption, Agung. The model sensitivity studies, however, do not support the 
contention by Ilirono et al. (41) that the discrepancy between measured and mo- 
deled Fuego decay rates can be accounted for by including aerosol growth in their 
model studies. A rapid growth process also is at odds with the Agung sulfate 
concentration data of Gan tinman, et: al. (43) which show maximum concentrations of 
sulfate occurring about one year after the eruption of the Agung volcano. Even 
if one assumed that the initial volcanic cloud contained no sulfate particles, 
one could only support the sulfate trends reported by Caatleman et al* (43) by 
assuming 1) a Large sulfuric, acid vapor to aerosol ratio at some time following 
the eruption and/or 2) that the bulk of the sulfate particulate mass is in a 
size range unaffected by sedimentation. For the first possibility to be true 
the sulfur dioxide to sulfuric acid gas phase chemistry would have to be much 
slower than estimated here. Both the chemical mechanism and the reaction rates 
in Appendix II could be incorrect. The second possibility would require a size 
distribution with fewer large particles for the sulfate portion of che aerosol 
layer and perhaps a reassessment of the adequacy of the present condensational 
growth mechanism. Measurements of the aerosol size distribution for a volcanic 
cloud in the stratosphere are sparse at best and are subject to largo uncertain- 
ties. It is also possible that the sampling of the stratospheric aerosol layer 
during the year after the Agung eruption was too poor to properly account for 
the expected layering and horizontal inhomogene ties and thus to establish with 
certainty the trends in total aerosol mass. 
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To summarize this section, H,,S0^ gas to aerosol conversion is show to 
occur rapidly. If the SO^ ’> gas phase conversion is essentially complete 

after about 10 da, then starting the growth model on October 29 is reasonable. 

As a result, the 1:1 gas to aerosol ratio estimated from Cadle et al. (2) for 
the volcanic cloud may be much smaller by the starting date (October 29) for our 
model runs. Therefore, the best estimate of the growth effects as shown in 
Figure 10 should be closer to the "no-growth" curve (circles) than to the solid 
curve. That best estimate is given by die dashed curve in Figure 10. At any 
rate, the 1/c decay time after January 1975 is relatively unchanged (2.5 months) 
from the case of no growth. This means that even though the particles gain mass 
by condensation, they grow larger and are subject to more rapid settling velocities 
and are therefore lost more quickly from the 16-21 km column. 
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5. DISCUSSION AND CONCLUSIONS 


Figure 13 describes the variation of the 16-21 km curves for 37N latitude. 

The lidar data in Figure 5 has been fit with a cubic spline function for easier 
comparison with the model result. The best estimate model aerosol curve is the 
dashed curve from Figure 10. The measured decay time t after January 1975 is 
8 months. The addition of each term in equation (l) yields: transport only — A. 5 

months; transport plus sedimentation— 2.5 months; and transport, sedimentation 
and growth — 2.5 months. The comparison in terms of Fuego aerosol layer decay 
times is not particularly good but it is believed to be within the uncertainties 
of the 2-dlmensional model calculation. 

By far the largest uncertainty (at least 100 percent) is due to the initial 
source function and the transport itself. This should not be surprising because 
it is difficult to obtain the zonal average of a point source function. Also the 
transport in the lower stratosphere at subtropical and tropical latitudes is 
particularly deficient because of the sparse data sets and inadequate theory. 

The sedimentation effect can alter t by up to 50 percent, because of the lack 
of good volcanic aerosol size distribution data and to a lesser extent the aerosol 
density. The sulfuric acid gas phase chemistry and aerosol growth terms do not 
seem to affect the results. This would not be true, however, if the gas phase 
chemsitry were much slower than the present 10 day estimate and if the 2.»rosol 
growth rates were to differ markedly from those for the background stratospheric 
aerosol layer. 

At the onset of this work, it was thought that a sensitivity analysis of 
the processes affecting the trends in the stratospheric aerosols was necessary 
to isolate the important deficiencies. That goal has been largly met. In ad- 
dition, it was felt that the Fuego event might yield a new set of "passive tracer" 
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BEST ESTIMATE" OF GROWTH 



Figure 13. Comparison of smoothed lidar measurements (16-21 km column number densities) 
over Hampton, VA with best estimate model results. 


data for actually validating the circulation in the lower stratosphere. This 
second goal has not been realised because of the uncertainties associated with 
the specification of the source function. In order to overcome this problem 
in the future, the model calculations can be initiated around January 1975 using 
all available aerosol measurements at that time and when the point source has 
been more evenly distributed around each latitude zone, 
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APPENDIX I 


Initial Size Distribution Estimate 

Stratospheric aerosol size distribution estimates have been critically 
reviewed by Harris and Rosen (34). Their data summary is applied to primarly 
non-volcanic material or to volcanic material that has decayed over time 
periods greater than one year. Five measurements of volcanic aerosols exist 
in the literature and they are 1) data for the time variation of the Agung 
dust from Mossiop (19); 2) estimates for the size of 0.15 to 0.25 micrometers 
from the photolectrlc particle counter of Hofmann and Rosen (8) ; 3) the data 
of Miranda and Dulchinos (35) for a January 1975 flight over New Mexico; 4) 
measurements of the St. Augustine volcanic plume by Hobbs, et al. (36); and 5) 
data on the volcanic plume of Mt. Baker by Rodke, et al. (37). Mossop used an 
itupactor measurement technique while all other data were obtained with some type 
of optical particle counter . Harris and Rosen (34) and Gras and Laby (38) have 
evaluated these two techniques with the conclusion that the in situ optical 
counter is more accurate and does not appreciably affect the physical state of 
the aerosol. 

The measurement by Hofmann and Rosen (8) on December 10, 1974 was processed 
to yield a Junge-type, dN/dr = Cr , (11) size distribution slope, v + 1, of 

3.1. This value is strictly applicable to the size range of 0.25 micrometers. 
That same slope was relatively unchanged on a second flight on February 18, 1975. 
Due to poor aerosol collection efficiency and the possibility that Mossop' s (19) 
Impactor may have been saturated with aerosol material from Agung, data from his 
impactor is not considered reliable, below a radius of 0.3 micrometers, so slopes 
are obtained from his data over only the range of -0.3 to 0.5 micrometer in radius. 
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His results vary from a slope of 1.0 on May 28, 1963 to 2.0 on August 6, 1963. 

On April 2, 1964, one year after the eruption of Agung, Mossop's slope parameter 
is 4.5. Data from St. Augustine have a mean value of 3.0 and a range of 2.6 to 
3.4. Data from Mt. Baker yield values closer to 2,0. 

Since most of the aerosol mass is concentrated in the larger particles with 
radii greater than 0.3 micrometers, it is that size range that must be estimated 
best. The eruption of Agung was at least a factor of ton stronger than Fuego, 
according to estimates by Cadle, et al. (2), while the eruptions of Augustine and 
Mt. Baker were even smaller than that of Fuego. Therefore, it is difficult to 
know what distribution might be most appropriate for Fuego. As a best estimate 
the August 6, 1963 particle size distribution of Mossop (19) was selected as 
representative of the initial condition for Fuego because its slope for particles 
greater than 0.3 micrometers in radius is 2.0 and that value matches the slope for 
the distribution in the Me. Baker plume. This initial distribution is plotted in 
Figure 3. Another distribution with a slope of 3.0 matches the estimate from data 
of Hofmann and Rosen obtained two months after the eruption of Fuego and it also 
agrees with the data of Hobbs, et al. (36) for St. Augustine, ?, less violent 
eruption. Clearly more measurement.* of aerosol size distributions in volcanic 
dust clouds are needed in the stages just following an eruption, in order to 
parameterize sedimentation and growth processes in model calculations. 
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APPENDIX II 


Gas Phase Chonistry 

In this appendix, a steady-state model of gas phase chemistry is described. 

The assumption made here is that during periods of non-volcanic activity the 

* 

background aerosol layer is maintained by conversions from gaseous sulfuric 
acid. 

To develop a chemical model for gaseous sulfuric acid, the following steps 
are necessary. 

1) Determine the possible reactions to be included in the chemical model; and, 

2) Determine background concentrations of all reacting species. 

A simple S0 2 to gas phase chemistry is assumed to follow the path 

so 2 -> hso 3 -> so 3 * h 2 so 4 

Background concentrations for these four gases must be calculated since observed 
profiles are unavailable. Since these reactions are all assumed to be one-way 
(i.e., H 2 S0 4 is the end-product of S0 2 reactions but H^SO^ does not convert back 
to S0 2 ) each conversion can be studied as a set of source and sink reactions. 

Park and London (22) have developed a photochemical model including reaction 
and photodissociation rates and species profiles. This particular chemical model 
was written as a subroutine to Louis' circulation model. The two models are com- 
bined in this work to do the chemical calculations, Louis' ’ circulation model is 
transformed into a one-dimensional model by by-passing the latitudinal computations 
to save computer time on these preliminary gas phase studies. As previously stated, 
the eddy diffusion parameters represent not only physical mixing in the atmosphere 
but also mathematical mixing that occurs when zonal averaging is performed. In 
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transforming the two-dimensional model, the eddy diffusion parameters reflect 

the additional mixing due to latitudinal averaging. The vertical diffusion co- 

3 5 2 -1 

efficients are assumed to vary from 10 - 10 cm sec 

The continuity equation for a trace gas is written as 
N(t) « N(t Q ) 4- [ ON/3t) ch + ON/8t) Cr ] At 
and solved for each time Increment. The approach is to assume an N(t Q ) and run 
the model until the quantity N(t) shows only seasonal variations. Then the con- 
tinuity equation reduces to (3N/3t) , * (3N/3t). indicating that the chemical 

cn tr 

and transport terms have reached equilibrium. The step At is chosen such that it 
is no greater than either the transport relaxation time (see Section 2) or the 
photochemical relaxation time. Since the chemical processes are faster than the 
transport processes, time steps are chosen such that they are about one tenth of 
the chemical relaxation time for the rate-limiting reaction. 

To obtain background concentrations of S0 2 , reactions (1) to (3) in Table 
II- 1 are considered. The chemical term is then written 
3S0 2 /3t = - k 1 (S0 2 ) (0H>’ (M) - k 2 (S0 2 ) (0) (M) 

- k 3 (so 2 ) (ho 2 ) • 

The dominant reaction can be determined by calculating the effective reaction rate 

^ A 

k according to.k.^ = k^ (OH) (M) . For the following molecular densities (22) 

OH = 1.5 x 10 6 cm ^ (at 20 km) 

and 

M = 1.8 x lO 1 ^ cm ^ (at 20 km) 

then 

k* = 2.2 x 10~ 6 sec" 1 (at 20 km, T ~ 230 K) • 
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Similarly, for values of 0 * 2 x 10 cm (22) and H0 2 x 10 cm”* 5 (22), k 2 * 

k 2 (0) (M) * 2.7 x 10 8 sec 1 and k!j " k 3 (H0 2 ) * 1.8 x IQ -8 sec” 3 . Reaction 

(1) is the rate determining reaction in this chemical scheme. The three e£~ 

it ik it 

fectlve reaction rates, k^, k 2 , and k^ are plotted versus altitude in Figure 
II-l. The net reaction rates can vary considerably with altitude either due 
to the change in the rate coefficients with temperature and/or to the change 
in the density of the reactant gases with altitude (see Figure 11-2) » Addl- 
tional errors can be introduced into the calculations of the net reaction 
rates due to the uncertainties in the rate coefficients and the diurnal change 
in density of the reactant gases. 

Figure II- 3 shows the resulting S0 2 profiles after running the model for 
100 days using the one-dimensional transport and photochemistry previously dis- 
cussed. The photochemical model of Park and London supplied the species profiles 
for the reacting gases OH, M, and H0 2> The S0 2 profiles converge to the 100 day 
profile as the photochemical and transport terms approach equilibrium. The model 

is run for another year to verify that any subsequent variations are seasonal. 

-3 

Jaeschke, et al. (10) have measured a value of S0 2 of 145 ng m STP for 13 km 
altitude. This value corresponds to a mixing ratio of 5 x 10 which compares 
very well with the 4 x 10 ^ mixing ratio at 13 km predicted by this model. 

The formation of HSO^ in this model occurs through reaction (1) . Possible 
destruction reactions are the following 

HS0 3 + OH + M H 2 S0 4 + M (4) 

HS0 3 + OH -> H 2 0 + S0 3 + M (5) 

Reaction rates have not been determined for the HS0 3 destruction terms; 
however, they are thought to be very rapid (29). Therefore, the net rate of 
change for the HS0 3 chemistry is assumed to be zero: 
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Figure H-1. Effective reaction rates of three SO destruction reactions. 




Figure II— 3. Model predictions of S0» for various days after start of rcodel 




3(HS0 3 )/3t - 0 - k x ($0 2 ) (OH) (M) - (HSO ) (OH) (M) 

~k 5 (HS0 3 ) (OH) (M) 

which after rearranging terms yields 

k 4 (HS0 3 ) (OH) (H) » k 1 (S0 2 ) (OH) (M) - k $ (HSOg) (OH) (M) 

The S0 3 phase of aerosol gas chemistry can be described by the S0 3 formation 
reactions (2), (3), and (5) and by the destruction reaction (6) from Table II-l 

SO^ + H 2 0 -> H 2 S0 4 k 6 « 9 x 10~ 13 cm 3 mol" 1 sec" 1 

The destruction term has an effective rate of 
» 

k* = k 6 (H 2 0) = 9 x 10" 13 x 10 -12 = .51 sec" 1 

and leads to a destruction of S0 3 that is fast compared to the effective S0 3 
formation rates, k 2 and k 3> As in the HS0 3 chemistry, 

3(S0 3 )/3t » 0 « k 2 (S0 2 ) (0) (M) + k 3 (S0 2 ) (HOg) 

+ k 5 (HS0 3 ) (OH) (M) - (S0 3 ) (H 2 0) 

or 

k 6 (so 3 ) (H 2 0) = k 2 (so 2 ) (0) (M) + k 3 (so 2 ) (ho 2 ) 

+ k 5 (HS0 3 ) (OH) (M) 

The final stage in this particular gas phase chemical scheme is the conver- 
sion to sulfuric acid vapor. Reactions (4) and (6) represent the formation terms, 
and an acid vapor loss rate was determined by imposing boundary conditions from 
vapor pressure arguments by Hamill, et al. (30) , 
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3(H 2 S0 4 ) - K 4 (hso 3 ) (OH) (OH) (M) + k 6 (S0 3 ) (H 2 0) 

The loss rate at the lower boundary, 8 km, of this one-dimensional model, id 

-12 

determined by setting (^SO^) equal to a mixing ratio of 6 x 10 , an equi- 

librium condition between condensed and vapor phase sulfuric acid. This 

causes a verticle concentration gradient at the lower boundary. The top of 

-12 

the profile, 50 km, was determined by setting Oh, SO.) equal to 1 x 10 
from previous equations the following substitutions can be made 

’3(H 2 S0 4 )/3t = k x (S0 2 ) (OH) (M) + k 2 (SOj,) (0) (H) 

+ k 3 (S0 2 ) (H0 2 ) 

Reactions (1), (2), and (3) determine the H 2 S0 4 profile shown in Figure II-4, 
with the exception of the sharp dropoff below 13 km and above 30 km, both of 
which are due to the imposed boundary conditions. A total (vapor + aerosol) 
sulfuric acid profile has been tabulated for comparison from data supplied by 
Hamill (27). The agreement is reasonable below 30 kilometers. That total 
curve in Figure II-4 is further broken down into vapor and aerosol profiles 
to point out that below 30 km all of the sulfuric acid is in condensed form. 

As can be seen from Figure II-l, the rate determining reaction for the 

H 2 S0 4 chemistry is reaction (1) (up to 40 km). The effective rate of this 

reaction is k^ * 2 x 10 sec which corresponds to a relaxation time of 
-5 

5 x 10 sec or about one week. In calculating this relaxation time, one 
must consider a factor of 2 error for the OH number density. Additional 
factors of 1.3 for the variation of the number density M r.nd of 2 for the 
variation in the rate k^ must be considered for a 10 km layer in the lower 
stratosphere. This means that the relaxation time for the S0 2 -> H 2 S0 4 con- 
version could be as much as 50 days. 
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H 2 S0 4 MIXING RATIO ppbv 


Figure II-4. Model predictions of H„SO, . 



The condensation of oulfuric acid vapor can be expressed by 

a 

h 2 so 4 + h 2 0 * h 2 so 4 • h 2 0 

and after n such steps 

H 2 S0 4 • (n-1) H 2 0 + H 2 0 -► H 2 SO a • nH 2 0 

Assuming n is 5 or less, the H 2 SQ 4 predictions, in vapor state, can be compared 
with aerosol observations (see the vapor and aerosol curves in Figure II-4 from 
data supplied by Hamill (27). 

Some features of this model which play an important role in determining 

the predicted H 2 S0 4 profile are presented in Table IT-2. The resulting H 2 S0 4 

-3 

peak concentration (and altitude of the peak) is given in molecules cm . The 
table also shows the assumed values and results used by three other similar 
steady-state 1-D chemical models. The conversion factors used to determine 
number densities from che various models are listed in the footnotes. 

The parameters presented in Table II-2 affect the resulting sulfuric acid 
profile in various ways. A larger vertical diffusion coefficient, K^, will re- 
duce the amount of sulfuric acid produced. The lower boundary condition, which 

affects the loss rate, and the K profile will have a major impact on the result- 

z 

ing altitude for the peak aerosol layer concentration. The dominant reaction 
S0 2 + OH + M will also affect the altitude of the peak due to the fact that the 
effective rate is a function of the S0 2> OH, and M profiles. Clearly, S0 2 is 
the source of sulfuric acid in these models and the amount of sulfuric acid 
produced will be directly proportional to the magnitude of the S0 2 concentration. 
The general results of this model then compare well with other published results. 
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Crufcsen (3) suggests that CSO photodissociation in the stratosphere contri- 
butes to the stratospheric S0 2 content. However, the CSO chemistry is only in- 
voked to explain the presence of S0 2 in the lower stratosphere in the absence of 
volcanic activity. It does not affect the ultimate conversion of SQ 2 to H 2 SO^. 

In Table II-3, the predicted values of F ,S0^ and SQ^ at 20 km are presented 

sat 

as well as observed concentrations of SO^ and particle mass concentrations. The 
most extensive observations are those of Lazrus and Gandrud (17) who have made 
near global measurements of sulfate mass mixing ratios. Their observed SO* can 
easily be compared to this model’s volume mixing ratios by multiplying by the 
ratio of the molecular weight of air to the molecular weight of SO^ (29/96). 

The 1971 values shown in the table are comparable to their more recent 1974 
background values. 


All values in Table II-3 have been converted to mass concentrations. The 

particle concentrations N , were converted to mass concentrations 

3 -3 

M by M =» N irpr 4/3 where valued of r * 0.3 micrometers and p » 1.5 g cm were 
P 

assumed. The volume mixing ratio in this model was converted to molecular number 
density and to mass densities according to M = (N/N ) x 98 where N is Avogadro’s 

a a 

number and 98 amu is the atomic mass unit weight of sulfuric acid. Aerosol con- 
centrations from Whitten’s model (29) have not been tabulated, since the percentage 
of sulfuric acid that is present in the aerosol phase is not given in his work. 


The discrepancy between the 1960 observations of Junge, et al. and of Friend 
and the later ones reported by Lazrus and Gandrud is not well understood. Lazrus 
and Gandrud suggest that one possibility may be that the impactor methods used in 
the earlier samplings are not as efficient in the collection of Aitken-sized par- 
ticles as the filter type method used in their own studies. Of course, another 
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possibility is that the differences are real. 


In conclusion, results of the steady-state chemistry for sulfuric acid 
or SO* fall within or close to the range measured by Laztus and Gandrud for 
sulfate at 20 km during times of low volcanic activity. An order of magnitude 
estimate of the relaxation time for the vapor phase chemistry, SOg l^SO^, 
then Is 10 days. From Table II-2, It is Important to note the various combi- 
nations of parameters which produce the results that are similar to the sulfate 
observations. Since measurement programs to confirm the proposed gas phase 
chemistry are lacking, it would be of Interest in future studies to determine 
the limits of cne various parameters (within the uncertainties) that would 
still result in reasonable comparisons with background sulfate observations. 
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Comparative Hodel Features and Corresponding Results 
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Table IX- 3 


Comparative Results Between Model Predictions 



Harrison and Larson (7) 
Harker (6) 

This model 


g cm- 3 
2.6 x 10" 13 

6.0 x 10~ U 


7.1 x 10 


-15 


i 


Measurements 


Junge, et al. , (1961) (11) 
(particle concentration) 

Junge and Manson (1961) (12) 
(SO* mass concentrations) 

Friend (1966) (5) 

(particle concentrations) 


Lazrus (1971) (16) 

(SO^ mass concentrations) 
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